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We report on the structural properties of Gei-xMux layers grown by molecular beam epitaxy 
In these layers, nanocolumns with a high Mn content are embedded in an almost-pure Ge matrix. 
We have used grazing-incidence X-ray scattering, atomic force and transmission electron microscopy 
to study the structural properties of the columns. We demonstrate how the elastic deformation of 
the matrix (as calculated using atomistic simulations) around the columns, as well as the average 
inter-column distance can account for the shape of the diffusion around Bragg peaks. 
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I. INTRODUCTION 

Ferromagnetic semiconductors have been extensively 
studied over the last decade as they are considered as the 
material solution for the needs of the spin-based electron- 
ics. Among the possible candidate systems, germanium- 
based ones have the advantage of being compatible with 
existing main-stream silicon technology. In particular, 
early results on the control of the ferromagnetism order- 
ing in gated germanium-manganese structures using an 
electric field^ have spurred many studies, focusing either 
on the diluted^ ^ or the heterogeneous^^ aspects of the 
GeMn system. Considering epitaxially grown GeMn sys- 
tems, several groups have observed the self-assembly of 
Mn-ri ch, columnar -like nano-objects embedded in a Ge 
matrij^l -*^^ * -*^^ ! ^^ * ^ -^. Curie temperatures above room tem- 
perature have been observed by Jamet et al. in these 
nanocolumns^^, and also more recently by Cho et al. in 
GeMn nanowireJ^, by Zeng et al. in homogeneous Mn- 
doped Ge thin films with xmu = 0.25%^ and by Xiu et 
al. in Mno.05Geo.95 quantum dot^, thus reinforcing the 
interest in this system. 

The low-temperature molecular beam epitaxy of Ge 
and Mn allows for a high concentration of Mn (up to 
10%) to be incorporated in the Ge matrix. Due to 
a spinodal decomposition occurring at the early stages 
of the growth, Mn-rich regions form. Depending on 
the growth parameters (substrate temperature, growth 
rate, manganese concentration), these inclusions can 
be spherical,!^ cigar-shapecP^ o r colum nar and extend- 
ing throughout the GeMn layei^^^E^. If the growth 
is performed at higher temperature or if the film is 
annealed, more stable phases can be seen, such as 
MnsGes, MusGes or MnnCeg clusters .^^lillllllllH^ 
Many studies have shed light on the magnetic and 
electronic properties of the GeMn nanocolumns and 
inclusions! ^^ ' ^^^^^ ' ^^ ' ^^ ' ^^ ' ^^^ ^ In terms of structural 



properties, the relationship between the GeMn la yer 
and the MusGea clusters is now well known. I^l ^^ * ^^ * Al- 
though information about the local environment around 
Mn atoms has already been reported, the understand- 
ing of the structural properties of continuous GeMn 
nanocolumns and their surrounding matrix is not com- 
plete. 

In this work we report on the strain and correlations 
of GeMn nanocolumns embedded in a Ge matrix, as de- 
termined by grazing-incidence X-ray scattering (GIXS) 
techniques at the European Synchrotron Radiation Fa- 
cility (ESRF), as well as elastic and atomistic simula- 
tions. While nanocolumns with diameters ranging from 
1.5 to 6 nm can be synthesized, this article focuses on 
samples with larger nanocolumns, where the inner part 
of the columns are either disordered or amorphous^^. 

The article is organized as follows: in part II we will 
present the experimental techniques used for the synthe- 
sis and the characterization (atomic force and transmis- 
sion electron microscopy as well as X-ray scattering); in 
part III we will show, using elastic and atomistic simu- 
lations, how the GIXS maps can be interpreted in terms 
of strain and correlations of the GeMn nanocolumns and 
the surrounding Ge matrix. 



II. EXPERIMENTS 

A. Sample preparation 

The (Ge,Mn) thin films were obtained using solid 
sources molecular beam epitaxy on epi-ready Ge(OOl) 
wafers. First, a 40 nm- thick Ge buffer was grown at 
250 °C after thermal desorption of the native oxide. Then 
Ge and Mn atoms were co-evaporated in order to get a 60 
to 80 nm-thick(Ge,Mn) layer with an overall Mn concen- 
tration ranging from ~ 6 % to ~ 10 %, depending on the 
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FIG. 1: (a) Plane-view TEM micrograph, (c) surface topog- 
raphy measured by AFM, and (b),(d) the central zones of 
the corresponding Fourier transforms (FT), each with an in- 
set showing the radial distribution of the FT amplitude. All 
color scales are linear, with the full z-scale for the AFM image 
corresponding to 2.6 nm. In each FT image a more intense 
ring can be observed, with correlation distances equal to ~ 10 
and ^ 40 nm, respectively in the TEM (inter-column average 
spacing) and the AFM (surface roughness) images. 

relative ratio of the Mn and Ge evaporation flux. The 
growth temperature was set to ~100 °C and the growth 
rate was ~ 0.2 A.s~'^. 

These growth parameters allow for a two-dimensional 
spinodal decomposition to take place at the early stages 
of the growth. Due to the layer-by- layer growth mode, 
the resulting (Ge,Mn) films are heterogeneous and made 
of columnar structures (nanocolumns) embedded in a 
matrix. In these conditions, nanocolumns have diameters 
ranging from 3 to 5 nm and densities between 10 000 and 
20 000 /im~^. The chemical analysis in a transmission 
electron microscope has shown that the composition in 
the nanocolumns is close to Ge2Mn while the surround- 
ing matrix is almost pure Ge (<1% Mn). 

It was also observed that samples not kept in an dry, 
oxygen- free environment would slowly oxidize over time, 
with a full oxidation occurring naturally after ~ 2 years 
of exposure. The effects of partial oxidation on the mag- 
netic properties of the GeMn nanocolumns have already 
been reported.!^ 

B. Transmission electron microscopy 

The samples were observed in cross-section and plane 
view by transmission electron microscopy (TEM) using a 




Qau,{nm-') 



FIG. 2: (color online) Bragg-GISAXS (see text for details) 
map measured by an angular scan around the in-plane (220) 
Bragg reflection, using a Vantec position-sensitive detector, at 
E=6.5 keV and ai = 0.4° (intensity plotted on a logarithmic 
scale, with 8 contour levels per order of magnitude). Qang and 
Qz are respectively parallel and perpendicular to the surface. 
Three finite distances effects can be seen, symmetrically about 
the Qang = axis : (i) a correlation distances of about 9 
to 15nm (marker 'D'), corresponding to the average distance 
between columns, (ii) the oscillations (marker 'T') along Qz 
(dot-dashed line) are due to the finite thickness of the sample 
(60 nm), and (iii) the two intense streaks (marker 'R') at very 
small Qang havc been attributed to the surface roughness. 

JEOL 4000EX microscope with an acceleration voltage of 
400 kV. Standard preparation, including mechanical pol- 
ishing and argon ion milling, was performed prior to the 
observations. Preparation for plane views also included 
wet etching using a H3PO4-H2O2 solution. 

The nanocolumns can be seen in plane-view TEM 
micrographs (Figjl^) as circular shaped objects with 
a different contrast from their surroundings. The 
nanocolumns induce strain in the matrix, as evidenced 
by the observatio n in high resolution TEM of a dark 
ring around them.^^^^ In some cases, the presence of 
a pair of dislocation^^ allowed for an estimation of the 
in-plane strain of about 4%. No such strain could be 
observed in cross-section TEM, where chemical contrast 
has shown that the columns are continuous, extending 
from the buffer layer to the surface. 

A fast-Fourier transform (FFT) analysis of the plane 
views reveals a higher intensity ring, as shown in 
Figjljb). Angular integration of the FFT as well as di- 
rect computation of the radial distribution function of 
the nanocolumns give an in-plane correlation distance 
between first neighbors of ^ 8 nm. 



C. Atomic Force Microscopy 

The surface of the samples was studied using atomic 
force microscopy (AFM). The surface presents a root- 
mean-square roughness trms = 0.8±0.3 nm, as mea- 
sured over 500 x 500 nm^ areas (Figflt) . A FFT analysis. 
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shown in FigjTJi, similar to that performed on the plane 
view TEM micrographs, exhibits a higher intensity ring 
at distances of about 40 nm, corresponding to a surface 
roughness correlation length. 



D. Grazing incidence X-ray scattering 

In order to get a better understanding of the correla- 
tions between the nanocolumns and their effects on the 
strain in the matrix, we performed X-ray scattering mea- 
surements on beamlines IDOl, BM32 and BM02 at the 
ESRF. We used energies around the Mn K edge (6.539 
keV) and the Ge K edge (11.103 keV). In order to en- 
hance the signal from the GeMn layer, we used a grazing 
incidence geometry, with the incident angle tuned around 
the total reflection angle {ac = 0.38° at the Mn K edge 
and 0.19° at the Ge K edge) in order to probe depths 
between 10 nm and 0.5 /im. 

The scattered intensity was measured using both Graz- 
ing Incidence Small Angle X-ray Scattering (GISAXS)^ 
as well as wide angle geometry (close to Bragg reflec- 
tions). GISAXS measurements are sensitive to the av- 
erage shape and positional correlations between objects 
at or close to the surface, but in our case contrast in the 
average electronic density between the columns and the 
surrounding matrix is very low, which prevents the obser- 
vation of the nanocolumn shape and distance correlation 
using this method. 

However a GISAXS-type analysis could still be per- 
formed around Bragg reflections: as the detector used 
for wide-angle measurements was a Vantec position- 
sensitive ID detector oriented perpendicularly to the 
sample plane, we measured the intensity as a function 
of the exit angle. In the case of an angular scan around a 
(220) reflection, as shown in Figj2j the 2D data presents 
features similar to classical GISAXS: the 'enhanced con- 
trast' is simply due to the sensitivity of GIXS near a 
Bragg reflection to any deformation (e.g. induced by the 
columns) in the layer. 

Two peaks can be seen in the 'Bragg-GISAXS' im- 
age: one (with d ^ 40 nm) corresponding to the surface 
roughness correlation (as shown in Figjlji, and already 
reported by Holy et al—) and another (with d ^ 8 nm) 
corresponds to the average distance between neighbor- 
ing columns, as shown in Fig{T]3. Moreover, the latter 
peak exhibits fringes with a period corresponding to the 
thickness of the GeMn layer (60 nm), which is only pos- 
sible because the columns are continuous over the whole 
thickness of the layer. 

In the large in-plane grazing-incidence X-ray scatter- 
ing maps, shown in Fig|3^, the most important feature 
is the absence of significant contributions far from the 
Bragg reflections of the standard diamond lattice. This 
is true for samples which are the focus of this article, 
with nanocolumns featuring a large (> 4 nm) average di- 
ameter, and either a disordered or an amorphous inner 
structure. ""^^ 
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FIG. 3: (color online) (a) Grazing-incidence reciprocal space 
map (RSM) of the measured intensity in the hkO plane - this 
map features almost only the (220) and (400) Ge Bragg peaks 
with a weak diffuse line between them, (b — e) high-resolution 
RSM around the (220), (400), (440) and (620) reflections. 
Around each Bragg peak an asymmetric ring can be seen at 
~ 0.07 reciprocal lattice units (r.l.u.): its positions is given by 
the average inter-column distance while the asymmetry of the 
intensity is due to the compression of the matrix surrounding 
each column. In (f-g) are shown the RSM around reflections 
(220) and (400), for the same type of sample as for Fig[2j 
but after oxidization (due to aging and exposure to air) of 
the sample: the correlation ring is almost circular due to the 
decrease of the strain of the Ge matrix around the oxidized 
columns. The intensity is represented as a logarithmic color 
scale, the full scale corresponding to 3 orders of magnitude 
for all plots. 
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A weak diffuse streak is also visible between the (220) 
and (040) reflections, with a faint diffuse (130) Bragg 
peak - these could either be due to diffuse scattering 
from the rough surface, or be related to the inner struc- 
ture of the nanocolumns (which is partially crystalline 
in the samples considered here) - this will not be the 
focus of this article, which mostly deals with the strain 
of the surrounding matrix and the correlations between 
nanocolumns. We will now focus on these characteristics, 
which can be analyzed using the shape of the diffusion 
around the germanium Bragg peaks. 



III. INTERPRETATION OF THE X-RAY 
SCATTERING MAPS NEAR BRAGG 
REFLECTIONS 

A specific diffusion is observed around all Bragg peaks, 
as can be seen in Fig 3]^b-e): this type of scattering 
pattern is due to the deformation of the diamond lat- 
tice, since the nanocolumns have a larger average lattice 
parameteil^^^ than that of bulk Ge. 

In FigjS] maps are shown for two different samples: 
Figjsja-e) correspond to the same sample with a strong 
asymmetry of the correlation ring, whereas Figjsjf-g) cor- 
respond to a similar sample measured after aging, and 
where the asymmetry is no longer present. This aging ef- 
fect has alreay been observed quantitatively using SQUID 
magnetometry,^, and has been linked to the sample ox- 
idation. In the case of a very old sample (^ 3 years in 
this case), it is expected that the columns are fully ox- 
idized. Although no direct observation (using TEM) of 
the nanocolumns morphology has been made, the strong 
decrease of the asymmetry can be related to the relax- 
ation of the matrix surrounding the nanocolumns, as we 
will now study. 

In the case of GeMn nanocolumns it was showri^ 
that the surrounding matrix is compressed. Therefore 
a Huang-type diffusioiP^^HH could be expected - how- 
ever the decrease of the intensity around a given re- 
flection {i.e. when following the intensity in a radial 
direction from the center of the Bragg position) is not 
monotonous, and features a maximum around a ring lo- 
cated at A{H or K) ^ 0.07 {r.l.u.)^ i.e. corresponding 
to a distance of ~ 8 nm in real space. This distance 
corresponds approximately to the average inter-column 
distance, and therefore the location of the ring is due to 
the inter-column distance correlation, while the polariza- 
tion of this ring (the intensity is larger on the high-q side, 
and also larger along [100] and [010] directions) is due to 
the deformation field of the matrix, as we will now show 
using elastic and atomistic simulations. 



A. Analytical model of the elastic deformations 

The scattering from the GeMn layer can be calcu- 
lated as the difference between the scattering of the dis- 




FIG. 4: (color online) (a) Amplitude of the atoms radial 
displacements (in nm), following an isotropic elastic model 
around a cylindrical inclusion, (b) Radial pair distribution 
function of the nanocolumns, measured experimentally (dots) 
and fitted using a gaussian + error function model, (c) Scat- 
tering from a single isolated nanocolumn around the (220) 
reflection, including the scattering contributions from the dis- 
placed atoms as well as the missing atoms inside the col- 
umn (see text for details), (d) Scattering of an assembly 
of nanocolumns around the (220) reflection, by taking into 
account their radial distribution function shown in (b). In- 
tensity is plotted with a logarithmic color scale for (c-d). 



torted and perfect {i.e. without deformation or columns) 
lattices^^- the advantage of this method is that it avoids 
most oscillation fringes due to the finite size of the simu- 
lated lattice, and also that it allows the superposition of 
the contributions from several insertions (columns in our 
case). The only approximation is that the sharp Bragg 
peak will not be obtained by this calculation. 

With the further approximation that the columns have 
a similar size and scatter independently, the scattering 
is equal to the interference of the scattering of all the 
columns: 



A{k) = Y,A,{k) 



,2i-Kk-Ri 



(1) 



where Ao(^) is the scattering due to a single column and 
the associated distortion of the lattice, k is the scattering 
vector and Ri is the position of column i in the sample. 
The intensity can then be written as: 

I{k) = \A{k)\^ = \Ao{mj2^''^'''^'\' 

i 

= \ Ao{k)\'^ X |7V + 2^cos27r^- (i?^- -i?^)| (2) 

i<j 

With the above approximations the scattered inten- 
sity should be proportional to the scattering of a single 
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nanocolumn, multiplied by an interference factor depend- 
ing on the pair distribution of the columns. If we assume 
that the inner part of the column is not contributing to 
the scattering (either disordered or amorphous), then the 
scattering of the single column (again calculated as a dif- 
ference with the scattering from a perfect infinite lattice) 
is the sum of the contribution from the displaced atoms 
in the Ge matrix and from the 'hole' left in the Ge matrix 
by the column: 

^o(^) = A{displaced atoms) + A{hole) 

displ.at. ^ hole 

= ^ fi e^'""^"^' ^e^i^k-u- - 1) - ^ f. 6^^^^-^° (3) 
* j 

where denotes the original position of atom and Ui 
its displacement. 

The so-called Huang scatterin^^^JHH typically yields 
nodal planes where the intensity is zero - which in this 
case would be expected for a plane perpendicular to the 
scattering vector. However this would only be true if only 
the A{displaced atoms) term was present - the contribu- 
tion from the hole makes the nodal plane disappear. 

We used the following analytical modeF^ for the radial 
displacements around a cylindrical nanocolumn: 

-(r > R) = + (4) 
47rr2(l - u) 

where R is the column's diameter, (5 is a parameter 
giving the extent of the deformation of the matrix, and 
u is the average Poisson coefficient for germanium. For 
the sake of simplicity we did not take into account the 
elastic anisotropy of germaniunP^ in this model. 

In the next section, we use atomic positions generated 
by this model to compute the scattering around Bragg 
reffections. 



B. Calculated X-ray scattering from the analytical 
model 

Grazing incidence X-ray scattering was calculated us- 
ing the atomic positions calculated from the elastic 
model: the displacements from an ideal germanium lat- 
tice are shown in Figj4^. To take into account the pair 
distribution of columns, we measured the pair distribu- 
tion using TEM plane view (Figjl^), and modeled it us- 
ing an isotropic radial distribution function, as can be 
seen in Figj4]3. As the sample exhibits no structure 
above the surface (embedded nanocolumns), transmis- 
sion coefficients^^ for the incoming and outgoing X-ray 
beams at the surface were not taken into account, as they 
only affect the scattering as a constant scale factor; the 
intensity was therefore directly calculated using Eq|2j 

The results of the simulation around a (220) reffection 
are shown in Figj4]3 for a single isolated column, and in 
Figj4]i when taking into account pair correlations. In 




FIG. 5: (color online) Atomic model obtained using a 
Stillinger- Weber potential, combining amorphous columns 
(grey spheres; darker shade indicates larger depth) sur- 
rounded by the germanium matrix (Ge atoms as black disks). 
Ge-Ge bonds of the perfect (undistorted) lattice are shown in 
green. Note the atomic displacements of the Ge atoms in the 
matrix, due to the column expansion. 

these simulations the intensity is larger at high Bragg 
angle, which is due to the compression of the Ge matrix 
by the nanocolumn. The pair correlation induces the 
correlation ring that is observed experimentally. 

However the experimentally observed intensity still dif- 
fers quantitatively from the simulated one, which is due 
to the ideal modeling of the lattice: (i) the columns are 
assumed to be identical (whereas the column diameter 
distribution is a Gaussian with an FWHM of almost 
~ 30%), (ii) the radial distribution will also differ de- 
pending on the column size, (iii) the true deformation 
of the Ge matrix should take into account the elastic 
anisotropy of germanium, and (iv) the distance be- 
tween columns (nearest neighbor distance of ~ 8 nm) 
is such that the displacement fields generated by all 
columns will not be independent. For this reason it is 
particularly interesting to perform atomistic simulations 
to obtain more realistic atomic positions taking into ac- 
count inter-column interactions, which is presented in the 
next section. 



C. Atomistic simulations 

In order to get a realistic elastic displacement field 
of the germanium matrix, atomistic simulations were 
conducted using a Stillinger- Weber potential^^ for the 
atomic interactions - one advantage of this potential be- 
ing that it was developed for both ordered and disordered 
condensed matter phases. To simulate a collection of sev- 
eral columns, we used a domain consisting of 100x100 
Ge-diamond unit cells, while the size along the column 
axis was 4 unit cells. Periodic boundary conditions were 



6 



used in the three space directions. The column density 
was chosen to be the experimental value, corresponding 
to an average of 42 columns per simulation cell (13 100 
/im~^). 

Starting from an empty simulation cell, the column's 
centers were added successively at random positions but 
with an exclusion distance between each new column and 
the previous ones. This exclusion distance had a Gaus- 
sian distribution with a mean equal to 6.9 nm and a stan- 
dard deviation equal to 0.35 nm. The resulting pair cor- 
relation has a peak at 7.6 nm close to the experimental 
value. The diameters of the columns were also set with 
a Gaussian distribution of mean and standard deviation 
respectively equal to 3.6 nm and 0.5 nm (corresponding 
to a FWHM of ~ 1.2 nm). Ge atoms were set on the 
germanium diamond lattice outside of the columns. 

Because we are primarily interested in this article by 
the matrix deformation and not by the precise structure 
of the GeMn columns themselves, we set atoms also inter- 
acting with the Stillinger- Weber potential, but at random 
positions inside the columns. Their initial density was 
the germanium density, but the lower density of the re- 
sulting amorphous structure leads to nanocolumns com- 
pressing the Ge matrix, as can be seen in FigjS] 

A standard minimization of the total potential en- 
ergy was then performed by conjugate gradient method. 
The resulting structure was a collection of amorphous 
nanocolumns embedded in a strained germanium lattice, 
as can be seen in Figs. |5]and[6^. 

To be able to simulate the incoherent diffraction from 
separate areas of a sample, we generated in the same way 
20 different domains (with the same statistical proper- 
ties) of 100 X 100 unit cells. 



D. Calculated X-ray scattering from the atomistic 
simulations 



Grazing incidence X-ray scattering was calculated us- 
ing the atomistic simulations: for the same reasons as be- 
fore, transmission coefficients for the incoming and out- 
going X-ray beams were not taken into account. The 
intensities were summed for all 20 independent domains 
of 100 X 100 unit cells, in order to average any 'speckle' 
structure that could arise from the absolute configuration 
of the 42 nanocolumns in each simulation. The result of 
these simulations are shown in Figj6]^b-c) for the (220) 
and (400) reflections, which compare well with the ob- 
served maps in Figj3]^b-c). 

The main difference with the results obtained using an 
elastic model (section III A) is that there are fewer cor- 



relation rings (only the first order is now visible), which 
is both due to the distribution of columns diameters in 
the atomistic simulation as well as atomic displacements 
which correctly take into account the influence of neigh- 
boring columns. 

In order to investigate the quantitative influence of the 
deformation of the Ge matrix on the asymmetry of the 
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FIG. 6: (color online) (a) Mapping of the simulated atomic 
displacements (from the reference perfect Ge lattice, displace- 
ments are given in nm), for a simulation with 100x100 unit 
cells and 42 GeMn nanocolumns (see text for details). Re- 
ciprocal space maps have been calculated by summing the in- 
tensities from 20 independent atomistic simulations, to lower 
speckle effects. The maps around (b) the (220) and (c) the 
(400) reflection exhibit the same features as in Fig|3jd,e). In 
(d) and (e) are presented the same maps as for (b) and (c), 
but where atomic displacements have been reduced by a factor 
5. The ring observed around the Bragg reflections is located 
at the position determined by the average distance between 
nanocolumns, and the polarization of its intensity (high vs low 
angle) is determined by the amplitude of the atomic displace- 
ments of the atoms in the matrix surrounding the columns. 
Intensity is plotted with a logarithmic color scale for (b-e), 
with a full scale corresponding to 4 orders of magnitude. 
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FIG. 7: Normalized angular distribution of the intensity on 
the ring surrounding the (220) reflection. This distribution is 
characteristic of the amplitude of deformation of the matrix 
lattice around the nanocolumns, with a larger asymmetry in 
the radial direction (0 vs 180°) for larger deformations: (a) 
distribution observed for a sample grown in situ (continuous 
line) and another (dashed line) with oxidized columns; (b) 
simulated distribution of intensities, using a single nanocol- 
umn, for various amplitudes of deformations - the relative 
displacement amplitudes around the columns are given in the 
inset - for the largest amplitude the maximum displacement 
is ~ 0.05 nm for atoms near the border of the columns. 



intensity in the correlation rings, we modified the sim- 
ulated atomic positions by decreasing the displacement 
(with respect to an undistorted Ge lattice) of all atoms by 
a factor between and 1. This approach is possible since 
atomic displacements are, in the elastic regime, propor- 
tional to the misfit introduced by the nanocolumns. The 
RSM maps for a factor equal to 0.2 is shown in FigjGjd- 
e), and reproduce well the rings observed for the sample 
with the oxidized columns in Figjsjf-g). 

A more quantitative comparison of the angular distri- 
bution of the intensity is presented in Figj7| the ampli- 
tude is maximum at ±45° with respect to the scattering 
vector, i.e. along the [100] and [010] directions around 
the (220) reflection, both in the simulation and in the 
experimental data. 



In FigjTjD, the angular distribution is shown as the am- 
plitude of the atomic displacements (with respect to the 
perfect Ge lattice) is reduced, from a maximum displace- 
ment of ~ 0.05 nm down to no displacement: there is a 
clear reduction of the asymmetry in the angular distribu- 
tion, which can be used as an indication of the amplitude 
of the Ge matrix strain around the columns. 

Note that the oxidizing mechanism occurring in the 
nanocolumns is not known here - obviously as the 
stress diminishes, the atomic density in the nanocolumns 
should also decrease, which probably indicates a migra- 
tion of some atoms of the columns to be replaced by oxy- 
gens. Also note that during this process, the diameter 
of the (oxidized) columns may have increased: but given 
the width of the size distribution of the nanocolumns 
(FWHM ^ 1.2 nm), this would only moderately affect 
the scattered amplitude. Therefore, the decrease in the 
angular asymmetry of the intensity can directly be linked 
to the decrease of the strain in the Ge matrix. 



IV. CONCLUSION AND SUMMARY 

In this work we have shown that grazing incidence X- 
ray scattering maps can be used for a quantitative anal- 
ysis of layers with GeMn nanocolumns embedded in ger- 
manium. In these samples the density of nanocolumns 
- related to the spinodal decomposition mechanism and 
the overall manganese concentration - leads to short 
(^ 8 nm) correlation distances, which can directly be 
measured from X-ray scattering maps under the form 
of correlation rings around Bragg peaks of the germa- 
nium matrix. Moreover the angular intensity distribution 
around these Bragg peaks shows a direct relation with 
the amplitude of the deformation in the matrix, which 
can range up to ~ 0.05 nm and is found to decrease 
upon oxidation of the columns. A more systematic study 
of the effect of aging in GeMn nanocolumns (including 
transmission electronic measurements and SQUID mag- 
netometry) is under way to fully understand the oxida- 
tion process and its effect on magnetic properties. 

Finally, while this article was focused on the study on 
the correlations between nanocolumns and the strain of 
the Ge matrix, several features indicate that there is a 
contribution for the inner part of the columns, such as 
the diffuse scattering between the (220) and (040) reflec- 
tions as well as the weak, large scattering around the 
forbidden (130) reflection, as can be seen in Figjsj This 
aspect of the structure of GeMn layers will be developed 
in a separate article, notably considering possible atomic 
structures within the nanocolumns l^^ ' ^^l 
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